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THE COROTATING VARIATION OP THE NORTH-SORTH 
ANISOTROPY OP COSMIC RAYS 

Xue Shunsheng Zhang Gongliang Xao Shaoyu 
(Institute of Space physics, Academia Sinica) 

Correlation analysis on the relation' of the 
north-south (N-S) anisotropy of cosmic rays, 
observed by the Nagoya multi-directional meson 
telescope, with the IMP as well as the solar 
wind velocity within solar Carrington rotation 
for the period 1971-1976. 

It is found that the N-S anisotropy of cosmic 
rays correlates quite well with the Bx component 
of the IMP. The correlation coefficient nearly 
equalsto 0.8. 

1. Introduction. Since the cosmic ray N-S anisotropy was 
discovered, a number of studies about it have been publidied. 
The drift due to the radial gradient of cosmic ray density n 
derives the N-S anisotropy in the presence of the IMP inten- 
sity B, as interpreted by SwinsontU . The correlation of the 
cosmic ray anisotropic index GG( [4] , [ 5 ]) with the daily 
average direction of the IMP was studied by Yasue et al. 1 ^ 1 , 
and a good correlation coefficient (~0.75) was obtained. But 
they did not give the quantitative analysis between the GG 
and magnetic field intensity. Fujimoto et al. l3J pointed out 
that the GG value decreased, as the IMP sector boundary 
changes its polarity from toward to away the sun, and vice 
verse. Munakata et al. C41 found that N-S anisotropy correlated 
somewhat with the solar wind velocity, as well as the Bx 
component of the IMF. 

We made the correlation analysis between the N-S anisot- 
ropy and solar flare, and found the flare of importance a 2 
with duration as 1 . 5 hours could change significantly the 
anisotropy 151. The present paper studies the relation of 
the corotating variation of cosmic ray N-S anisotropy with 
solar wind velocity, as well as the parameters of the IMP 
with the sun. 


2. Data and Method. The cosmic ray data used are the daily 
average values taken from the multi-directional meson teles- 
cope at Nagoya for the period 1971-197616]. The cosmic ray 
N-S anisotropic indexes A and AA are defined below, 

A® (N-S)/2 
AA=(N-S)+(N-E) 


a] 


Another parameter I is defined as 

I*(N+S)/2 (3) 

Each term from formula (l) to (3) denotes a directional 
component of the cosmic ray intensity having its central 
direction of the viewing cone pointing toward the zenith 
angle of 30° in the north (N), the south (s), and the east 
(E)-directions, different from that adopted by Munakata et 
al.. The solar wind velocity and the IMF data are taken from 
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that compiled by King, J. H. C7J . 

Analysis of cosmic ray data are made by applying Chree 
superposed epoch method with reference to their ordinal day 
in a solar rotation® . An accurate Carrimgton rotation 
period is divided into 27 round days. It avoids not only the 
drifting error caused by adopting Bartel’s rotation, but 
also the difficulty produced by decimal Carrington rotation. 
Such a converted round day is called the "Carrington Day", 
with length approximately equal to 1.01 calender days. After 
taken the nearest whole number, the day at which the zero 
Carrington longitude crossing the solar central meridian is 
regard as the beginning day of each Carrington rotation. 

In all figures of this paper the absissas represent the 
Carrington ordinal day from the first day to the 27th day 
corresponding the solar longitude from 360° to 0° . The 
ordinates show the parameters concerned, i.e. A, AA and I, 
as well as the solar wind velocity V, IMF total intensity 
B, its components Bx and By and its azimuthal angle in the 
ecliptic plane. In order to making the figures more clear, 
the values of the first three days are drawn repeatedly 
after the 27th day 


3. Analysis and results. 
N-S anisotropy and IMI 1 


The distributions of 
Bx component in the 
rotation have been calculated. Simply 



Ftg. I. Corototifl? vortotlons of (, A, AA ond Bx. 


cosmic ray 
Carrington 
shown are the results 
for 1973 and 1974 in 
Fig. 1. In Fig. 1b, 
two steady corotating 
structures of the IMF 
can obviously be seen* 
Also noted are the 
variations of A and 
AA being similar to 
that of Bx component. 
In other years, the 
steady corotating 
structures still exi- 
st, but not as clear 
as in 1974. Because 
the correlation coef- 
ficient between the 
A and AA is almost 
closed to 1 , only the 
AA is discussed in 
following. 

between AA and Bx and By 


The single and bi-correlation 
can be expressed as the following formulas, 

AA=a 0 , +a, * B x ^41 

AA=ao2+aiyBy (5) 

AA=ao+a 1 1 Bx+aiaBy (6) 

Table 1 summarizes the yearly coefficients in formula (4) to 
(6) and respective correlation coefficients rx» ry and r X y, 
as well as the I and AA. It is striking that AA correlates 
quite well with the Bx component with a high correlation 
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coefficient rx 0.8 The ry is negative, its magnitude is a 
bit less than rx* The bi-correlation coefficient rxy is a 
little higher than rx* The rx obtained here is much higher 
than that by Munakata et al.. This is because the present 
correlation coefficients are calculated on a yearly beses, 
while the previous ones using five year data altogether, 
and the change of the AA following the solar cycle reduced 


Table 1 . Values of various coefficients, as well 
as the I and AA during 1971-1975 


year 

a oi 

a 

B 

Q 

a 02 


Efl 

ao 

a 

1 1 

a 

12 

r X y 

I 

AA 

1971 

14 

4 

.2 

0 . 

.86 

14 

-4.1 

-0, 

.80 

14 

3 

.0 

-1 

.4 

0 

.87 

1705 

3 

1972 

38 

4 

.2 

0 . 

.71 

39 

-4.1 

- 0 , 

.63 

39 

3 

.1 

-1 

.7 

0 

.74 

1557 

38 

1973 

69 

4 

.0 

0 , 

• 73 

70 

-4.3 

- 0 , 

.69 

69 

3 

.4 

-0 

.8 

0 

.73 

1400 

64 

1974 

90 

3 

.4 

0 . 

.88 

92 

-3.6 

-0, 

.88 

91 

1 

.8 

-1 

.6 

0 

.88 

1294 

91 

1975 

114 

3 

.3 

0, 

.82 

113 

-3.3 

-0, 

.78 

114 

3 

.5 

-0 

.3 

0 

.82 

1187 

97 


the statistical accuracy. Meanwhile the superposed epoch 
method has filtered the short term vaviation and raised also 
the statistical accuracy. The coefficients a>x , aiy , an and 
a t 2 do not change very much, but aoi , ao 2 and ao increase 
approximately from 14 to 114 with the decreasing solar 
activity during 1971-1975. This fact illustrates that the 
variation of N-S anisotropy is due to not only Bxvn, but 

other factor, for example, the variation soTar general 
magnetic field. 

We have also calculated the distributions of cosmic ray 
N-S anisotropy, solar wind velocity and the IMF parameters 

in the Carrington 
rotation for this pe- 
riod. For simplicity 
_ “ T only the results for 
0 1971 and 1974 are 
~ 2 shown in ELg. 2. Table 
o , 2 summarizes the cor- 
« relation coefficients 
■ 6 of the AA with the V, 
¥ and B respectively. 
The magnitude of r v , 
““'•ry and r B are much 
lower than those of 
-500 Bx and By. Hence it 
is the best to study 
the correlation of AA 
with IMF component . 

Fig. 2. Corolating variations of A A, Bx, By, B,y ond V. 

It is well known that the solar activity weekens from 

1971-1974, the solar modulation of cosmic rays increases, 
and the omnidirectional intensity of cosmic rays increases 
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Table 2. The ry» Ty and 
r B for 1971 and 1974 


Year 

rv 

Ty 

r B 

1971 

-0.33 

-0.61 

0.29 

1974 

0.21 

-0.75 

0.02 
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every year. The present results 
show the AA follows a similar 
trend, however, it is unexpected 
that the I decreases from year 
to year, quite different from 
the variation of omnidirectional 
intensity. 


4» Conclusion. The cosmic ray N-S anisotropy correlates 
quite well with the IMP components Bx and By with a correla- 
tion coefficient nearly equal to 0.8, a positive one for 
the Bx component and negative for the By. The bi-correlation 
coefficient with both Bx and By is almost closed to 0.8. The 
variation of the regressive coefficients illustrates that 
the cosmic ray N-S anisotropy, in addition to the gradient 
drift effect, is controlled by some kinds of long term 
varlayion of the solar cycle. 
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